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I 

Oh iht rapid calculation of. times of moon-rise and moon-set. By J. H. Field, M.A., 
Imperial Meteorologist, and S. M. JACOB, I.C.S. 

There has arisen of late years in India a considerable demand on the Meteorological 
Department for information regarding the times of rise and set of sun and moon. 
During, a large part of each year the weather conditions in this country are sufficiently 
•definite to allow of the construction of street lighting tables, based on sunlight and' 
moonlight, and it has accordingly come to be realised in large towns that a great saving 
of oil or other illuminant can be arranged for when precise information is available as to 
the hours of darkness to be expected. 

3. As long as it was required to make the calculations for one place only, the 
routine method described in paragraphs 3 to 6 was adopted.* When, however, the neces- 
sity arose of computing times of set and rise of the moon for many places throughout each 
year, the work involved became sufficiently great to call for shortened methods ; and the 
principal object of this note is to show (paragraphs 7 to 1 1) how, when data for the moon) 
have once been calculated for a station A by an unskilled computer, the calculation can 
be subsequently, extended to any other geographical positions B, C, etc., with compara- 
tively little further expenditure of time. 

.•3.. For the routine determinations for station A it is first necessary to calculate 
a- table of hour-angles in terms of declination (1. e., of time intervals from 
meridian pa.ssage for all declinations), and then to apply the particular values of hour- 
angle, relating to the declination at each rise and set, to the daily times of meridian 
transit, separately calculated. 

• This preliminary routine method, paragraphs 3 to 6, probably contains nothing intrinsically novel, bat has been • 
described bcc.ansc It hfls proved satisfactory where only unskilled computation is aTallable, subject merely to ucrotiny of 
nnlK for consistency, and to examination of the working for a day here and there. 
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The hour-angle declination tabic is obtained as follows 

In t P is the pole, i the xenilli, Q 0 the equator, b h the 
horiion, m the moonl 

Then Pht=latitudc 

o =declination, + when north. 

2= zenith distance. 

tlsrhour-angle. . 

Pm=ga — o. 

Pz=go — 



Therefore in the triangle Pm z, cos Z 


cos 


g — cos Z— sin 5 sin 


sin S sin ^ + cos S cos ^ cos 6 
i:' 


COS B cos 6 


The calculation of 0 for a given place results, for the latitude of Allahabad, in 
a curve such as fig. 2 : Z has there been taken as 90° minus 45' to allow for refrac- 
tion, parallax, and for a position of the moon such that its lower limb just touches the 
horizon.® This is the most convenient position to consider, since it is then unnecessary 
to go into the question of the orientation of the illuminated sector, which varies with 
the moon's age. 


4. With the hour-angle table ready, it now only remains to caltulale the daily times 
of meridian passage in local mean time, apply the hour angle proper to the declination 
for the time of set or rise, and then transform into civil time for the standard meridian 
of the country. 

Now the L. M.t time of meridian passage at any place differs from the G. M. T. of 
Grccnwicl) meridian passage merely by the amount of change of moon's R. A. in the 
time interval between the two meridians. The daily figures for Greenwich are given in 
N. A, (table IV for each month) and from this dally difference, due to daily change 
in R. A., there follows the amount of the correction required. 

It is necessary next to find from the curve, fig. a, the hour-angles to be applied 
to the time of meridian passage just found, and to do this the moon’s declination at 
rise and set must be known. The declination is given in the N. A. for each hour of each 
day (tables V to XII monthly) but appears there inG. M. A. ,T. It is' necessary 
therefore to determine approximately the limes of Iccxil rise and set, put these into 
G. M. A-T. .and then enter the declination from the rdnir.nac. For the present purposes 
tl-.e times need only be rij'-proximatc .since declination change.? slowly, and they, may be 
taken as 6 h.ours before .and after meridian p.a.ssage.f 
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• The sums of the time of meridian passage and the hour-angles have finally to be 
corrected for change of the moon’s R. A. in the latter intervals, and this is done by refer- 
ence to the N. A. for rate of change at the time, 

5. The actual process of computation is most conveniently done by means of the 
tabular statement below, in which the columns are to be filled in by the computer with 
the figures indicated in the explanatory notes below it. 



The example given is for Allahabad in long. Si° 54' E. lat., 25° 2S' N ; and the 
following notes apply : — 

Column I ; date. 

„ 2; Greenwich Mean Asiron. Time of Greenwich upper meridian passage, to 

nearest minute, direct from table IV monthly in nautical almanac. 

„ 3 ; Change of column 2 in 48 hours, to nearest minute, 

„ 4; Col. 3x0/4$ where n= hours ca'l of Greenwich: t'n4 for Allahabad). 

, „ 5 ; Col. 2 — col. 4; local mean asltoii. time of loc.il meridian p,ass.ige. 

„ 6 ; G. M. A. T. of local moon set or rise (sufTiclcntly accurate if put equal to col. 5 

in G. M. A. T. +and — 6 hours ;-}-for set nnd— for rise). 

,, 7; J.foon’.s declination (nearest i”) at lime of col. 6 , 

„ S; Hour-angle (to i min.), from curve of 0 , S; phi-s sign for set, minus sign 

for rise. 

„ 9; Correction (to ( min.), for change in R. A. of moan during time col. S round 

abo:it time of col. 5. This latter should really be put into G. M. A. T., but 
a.s R. A. ch.inges slowly, the error in using col. 5 is not more than 4 min. 
for any longitude. Signs .as in coi. S. 

„ to; Col. S -I- Col. 9; lapse of time .since meridian passage. 

,, ti ; Col. 54- Col. to : L. M A. T. of local rise or set. 

„ 12 ; Col. It + 12 hour.s, to reduce to Civil Time, + correction (for Allaliabad =- 1 - 3 ^ 

minutes) to the nearest sl.indard meridian, S2^°L. for India. 

> correction of dale, :.c., 27 li. 55’3 m. on 9th becomes 3 li. 55‘2 m. on totli. 

>. Mi) 
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6, With the approximations employed in the tabular form, the error in the rarious 
steps will not 'exceed, for India, . 


Column a 

4 »l 


± 

o'S 

minute 

Columns 3 and 4; 

• «« 


■± , 

0*1 

tr 

Column 6- 



± 

I'O 

n 

u 7 

4t« 

• • • 

± 

0*7 

fi 

s 




0*35 

i» 

.. 9 

i«4 

•ft 

± 

o*as 

ej. 

f} - 


or the maximum possible accumulated error will be ± 3 * 8 ; min.' in the; times of rise: 
and setr 

Method of cxtcQjdinff the results of Station A to a number of 
otherlgcographical positions. 

•j. The additional process now to be described, is used 'for all further places ■. 
requiring moonlight information, and is based on the fact that it is easier to transform • 
times from a station, A, preferably on the equator® to a series of other geographical posi-- 
tions B, C, etc., than to calculate each of the latter separately by the process already t 
described. 


With the notation employed in paragraph 3 we have 
cos Z 


cos 


e = 


sin <i sin S 
^ = a, say 


cos 0 cos S 


so that 

0 , — cos —'a — cos b ~ cos~'k, say 

k = cos (cos -'n — cos “'i) 

= ab + sin (cos ~'a) sin (co 5 “'i) 

= -f v/“(i — «’) (1 — 6’) 

0 , -~e, = cos - [ fli + v/~( 1 _ a ') ] 

When Z = 90°, as in the case of a star, with refraction left out of account, we 
have : — 

6 , 0 , — cos“‘ [ tan 4 >, tan •}>, tan’S 

+ (i — ^^tan’ 4 >, tan ’S) (1 — tan’ <f), tan 'Sf] 
put = o, i.c, assume a place on the equator, 

then 6 , — 0 , = cos~’a./ t — tan’ </», tan ’S 
= sin (tan tan 8). 
sin ( 0 , — ) =:tan 6 . tnn 8 

a result c.ailir.g for the application of a small correction, paragraph 10, when Z is not 
made exactly 90^. 

8. This equation has been worked out for various latitudes, each in steps of 5' of 
declination, giving 8, >— 0 , in inierv.als of sidereal time. After correction to mean-time 
intcn'als, and allowing for the change in R. A. of the moon during those intervals 
(at an tivet.ige rate of s's sees, in 1 min., an average which gives a maximum error of 
I min. for latitudes less than 35®), the following figures for d, — 0 , are obtained : — 


• YlS =-'* sn j'lt-c-n Iff ll e the tr.it to the aprrclir.a'.isn the 

ra i-i'k:jcci;n sve rcee -whleh foi tree.-. rraJiM eoksias! .-.t 5 heart 57 cirstet. stid jer the 

»',-e «rr»ft et 6 1: -a-i ; njitet-r. risr iriJ jet leirp O-Cfed ss ir. paiapsph 3 ss-; ;*• loelr.o'c. 

Y!{-?c»rr<t.-tT-y rir-pVTicnti-n Hthc pr.mrr.'.pb 5. tiecc llie v.-iri.'.litc c! cotorr.n 8 

fcecnrrtj RCc-ftT',, t-'or tH tun tre jTrrS,f;£.-.-.;nr, h -ti'l prc.Wtr,.nfd iV'c'cn'ciifntion cf Tire ‘ent! sit in L. M.T, (or 
ery f -r'i'."' nervit cnrOjtj ir .'.npifiri; the e(juo!i{:r. c; erre.nr.d itr rcrrcctior.r. fCT l.ititudc (fre-r. figures) to the 
cc,-'! -nt K"';."..nr;:!e f I inn; 3 mir.u’.ri. 
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These nrc shown plotted in fig. 3. 

9. The trnnsfonnation is most conveniently made wlien .an equatorial position has 
been chosen for the place of reference A, but if for any re.ison the original figures relate 
to another latitude Allahabad in the example above) the transformation can 

obviously be cfTccied by employing the di^crcucc between the corrections, from fig. 3^ 
towards the equator and away from it to the new situations under treatment. An illus* 
tration of the application of the transformation table may bo useful. Taking Allahabad, 
lat. 35° 28'N long. 81“ 54'E, as the sMtion of reference, and starting with col. 1 a as 
previously determined in paragraph 5, the times for Pcslnwar, lat. 34* I'N, long. 71° 
SS'E, m.ay be obtained as shown in the following table: — 



Ctl. 7 



SO 

ss 

S3 


A 

n 

C 

D 

ul 

pih 

li'N 




1 

3r-55'3 


+ 1 1*5 

28-497 

4-497 

loth 

riS0 






rr 

i —14 

! 

'5-3 7 


gth 


Here the columns 7 and ta in italics .arc those already worked out, i« e.v/cHro for 
Allahabad, and the blank table for the new .station is merely to he clipped to the table 
of figures for Allahabad while the remaining work, columns A to D, is done. Column 
A contains the latitude coTrccUoits from the curves, figure 3, for the declinations of 
column 7, and in the present case is made up of the difference of correction from 

Allahabad towards the equator, and away from it to Peshawar. 

' \ 

7 hus in set, declination N, wc imve — a8'5+40=-f ii’s min. 

.. rise. N, „ +35— -49=— *4 

(for declination S, the signs would he reversed.) 
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6 . With the approximations’ cmployed:in the tabular, form, thfi error in the verious 
steps ■will not ‘exceed, .for India,, • • • . / . ■ 


Column 3 

«•! 


±' . 

o’5 

minute 

CoIumns!3.and 4: 



db’ 

0*4 . 

11 

Column 61 


• «ft ' 

± 

I’O 


7 



.'±‘. 

o‘7' 

it ’ 

„ 8 


ftftft 


o'ss 

^ it . 

.. 9 ., 

• •ft ■ 

««• 
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o'3S 



or the maximum possible accumulated' error will be db 3*8 min. in the , times of rise ; 
and set. , ’ • , . , 

Method of extending the resiilts of Station A to unimxber of 
other\geograiyhica,l positiojis, . : 

7. The additional process now to be described, is used '.for all further places ^ 
requiring moonlight information, and is based on the fact that it is caster to transform * 
times from a station, A, preferably on the equator* to a scries of other geographical post-' 
tions B, C, etc., than to calculate each of the latter separately by the process already! 
described, ' - 

With the notation employed in paragraph 3 wc have 


cos 0 = 


cos Z — sin </> sin S 
cos 0 cos S 


a, say 


so that 


d, — = cos “'c — cos i = cos~’,t, say 

k = cos (cos ~*rt — cos ~'b) 

= fli + sin (cos ~'fl) sin (cos 

= flj -h \/~{J — a') (i — 6’) 

0 ,- 0 , — cos “' [ rti -i- 1 

When Z = 90°, as in the case of a star, with refraction left out of account, we. 
have : — ■ . 

6 , — 0 t~ cos"" [ tan tan (f>, tan’S ' . 

+ a/ (t — .tan’^, tan ’ 5 ) ( i — ’ lan' ij>, tan ’S)~) 

put = o, t.c. assume a place on the equator, 

then 0 , — 0 , =■ cos'”,/ i — tan’ ^i'tan ’S 

= sin (tan tan S). 

.V sin ( 0 , — t?, ) =tan t/i. tan S 

a result calling for the application of a snudl correction, ,, paragraph to, when Z is; not 
made exactly 90“, . 

8. This equation has been worked out for various ' latitudes, each in steps of '5' of 
declination, giving 6 , ■ — 0 , in intervals of sidereal time; After correction to rhe.an-lime 
intervals, and allowing for the change in R, A, of the moon .during those intervals 
(at an average .rate of a’e .sees. In- 1 min., an average which gives a maximum error of 
I min, for latitudes less than 35°), the following figures for 0 , — 0 , arc obtained : — 


• The use of an equatnrinl slnlion for the otinintl cntcn'.ation has the ndvanlano that lo the rerjuirr itanj'fo*''’'*'''’'’ 
declination has no itiHuenceon the hour ancle, rvhicli for the nioon ttnuins constant nt 5 hour.* 57 minutes, and for the 
tun constant at 6 hours 3 minutes, rise nod set being denned os in paragraph 3 and il, footnote. 

TWs constancy introduces n simplification in the 'c.n!cttlation 5 of parngrnpli 3, since liic v.iri.al)Sc of column 8 
become.? n constant. I'or the sun tlie 6inipIific.ntioti is Mill grc.itcr, and the ’c.nlc.ul,nticn cl risc'ond stt in I,. M. T. for 
any position merely consists in .applying the cqu.ation of time .and the corrections for . I.ilitutic (from ’figure 3) to the 
constant hour-angle 6 hours 3 minutes. 
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These are shown plotted in fig. 3. 


9. The transformation is most conveniently made when an equatorial position has 
been chosen for the place of reference A, but if for any reason the original figures relate 
to another latitude (c.g., Allahabad in the example above) the transformation can 
obviously be effected by employing the difference between the corrections, from fig. 3^ 
towards the equator and away from it to the new situations under treatment. An illus- 
tration of the application of the transformation table may be useful. Taking Allahabad, 
lat. 35° 28'N long. 81® ■54'E, as the station of reference, and starting with col. i a as 
previously determined in paragraph 5, the times for Peshawar, lat. 34* I'N, long. 71° 
3S'E, may be obtained as shown in the following table : — 
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t4-34'7 


—14 
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Here the columns 7 and i a in italics are those already worked out, tn exteuso for 
Allahabad, and the blank table for the new station is merely to be clipped to the table 
of figures for Allahabad while the remaining work, columns A to D, is done. Column 
A contains the latitude corrections from the curves, figure 3, for the declinations of 
column 7, and in the present case is made up of the difference of correction from 
Allahabad towards the equator, and away from it to Peshawar. 

Thus in set, declination N, we have — 38’5-|-4o=-fn’s min. 


„ rise, „ N, „ >f3S— 49=— 14 min. 

(for declination S, the signs would be reversed.) 
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Column B = column lo -}- column A + a correction, constant ihroughout the 
ycar^ composed as follows ■ 

1 

(rt) +4rs min. ithe difference between the meridians at Peshawar and 
Allahabad. 

(b) 4-r5 min. the correction for change of the moon’s R. A. during the 
interval (a). 

Columns C and D allow, ■ as before, for any change of d.ate resulting from times in 
column B c.\cccding 24 hours, and only parts of ihem would be filled in., 

10. Errors . — In the process of transformation as illustrated above, three small 
errors are added to those detailed for station A in paragraph 6. 

We have first the error due to the assumption that over the range of latitudes 
considered, the obliquity of rise and set has acon.sfant value, conveniently made equal 
to the mean value for the country in question. Second, the error arising from the change 
of declination of the moon during the interval 0 ,- 0 , for dificrcnccs of latitude, and 
during the time interval involved by dificrcnccs of longitude. Third, the error due to 
assuming in column B item ( 5 ) a uniform rate of change in R. .•\. of 3‘2 secs, per min. 
of time. 

To consider the first of these. The transformation curves have assumed Z=go*’, 
and if, as is convenient, it is intended to give times for ZreSg'ts', with the moon’s 
lower limb just toucliing the liorixon, a correction must be made, which varies to a 
slight e-xtent with latitude, but for practical purposes can be tahen .as constant over & 
wide range of latitude. 

If t be the minutes of lime for the moon, with declination 8, to rise a minutes of 
arc in latitude </> we liavc as usual 


iS\/cos’S — sin’'f' 

and the greatest and least values of this, when n=.)5', are 3’5 min. and 3*4 min. 
respectively, with the limits of latitude in India ; and 3 ‘3, 3'o the corresponding values 
on the equator. If then we have calculated tlic hour-angles for station with Z =r 
90° — 45', or if we apply once for all to the v.aliies calculated for Z=£;o° a correction of 
3 rain, throughout in column n, the neglected variation of this correction with latitude 
over the wide range required (say 34‘’N to 34'S), will introduce an additional error of 
only ±o'6 min. of time. 

In regard to the second error, it is to he noted that for a given place the change 
of declination in the time interval Oi-O, is allcrn.atcly plus .and minus for both positive 
and negative declination, while during the sm.aller interval in longitude the change in 
declination changes from plus to minus once only in each complete swing in declination. 
The variation, with moon’s declination, of this combined error is not therefore a siniple 
curve, and it is best in cases where I'S min. isaiot of importance in the result to con- 
sider it an error, and not to attempt to correct it. 
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The third error never exceeds o'8 min. for the whole range of longitude in India, 
when the equatorial place of reference has been chosen for a mean longitude. 

The maximum accumulated error cannot then exceed ± (2'8-}-ro+o‘6-l-i-5+'8) 
67 min., and will in general he very much less. 

II. With regard to the time taken for computation, it has been found that the 
complete year’s figures of the station of reference A (on the equator or elsewhere), can 
be worked out and checked in 1 8 hours (3 days) of a clerk's time, and that the sub- 
sequent transformations for stations B, C, etc., in other geographical positions in India 
need take only 6 hours each per year. It is in general unnecessary to deal with the 
figures for rise and set which occur during day-time, or even within i hour of day-light, 
and in some case the periods covering the last and first quarters can be disregarded ; the 
omission of these periods still further shortens the work. 



Diagrams. 

Figure 2 gives the hour'angle declination curve of the moon, when its lower limb 
just touches the horizon at rise and set, refraction and parallaK being allowed for. 

If the scales of time and declination arc reversed in direction, the curve gives the 
hour*angIe for the sun at first and last ray, refraction included. 


Figure 3 gives the correction necessary to obtain the hour*anglcs of rise and set in 
latitude >fy, from the corresponding hour-angles on the same meridian at the equator. 
The correction is positive when hour-angle and declination arc of similar sign, and 
negative otherwise. 
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